A New Indicator of Nonlinear Gravitational Clustering 
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Abstract. Alignment of velocity and acceleration before shell crossing, 
and later misalignment are used to define velocity contrast, an indicator 
of dynamical state of matter undergoing gravitational collapse. We use 
this to study bias in clustering properties of dynamically nonlinear mass. 



1. Introduction 

The problem of gravitational dynamics in an expanding universe can be re- 
duced to the corresponding Newtonian problem for a nonrelativistic coUisionless 
medium, if scales of interest are small compared to the bubble radius. Further 
if = ^Imatter = 1 then following equations describe the system. 



Here u = dx/da is peculiar velocity, g is "acceleration" and tp is gravitational 
potential. These quantities have been defined in a dynamically relevant manner, 
their relation with the usual definitions are given in (1). 

All concepts introduced here can be defined for a density and velocity field 
or its particle realization. In this paper, we use particle picture for convenience. 

Dynamical evolution described by (1) can be divided into four regimes. To 
motivate the form of new indicator, we describe evolution of trajectories in these 
regimes. In linear regime density contrast grows in proportion to scale factor. 
Peculiar velocities are constant and equal the acceleration, u(x, a) = g(x,a). 

This equality can be extended, approximately, beyond linear regime. Tra- 
jectories closely follow those predicted by Zeldovich approximation ( Zeldovich, 
1970) upto shell crossing. Particles move with constant peculiar velocity, veloc- 
ity being fixed at the initial position : u(x, a) = g(xj„, aj„). Accuracy of this 
approximation implies u(x, a) ~ g(x, a) for the Zeldovich regime. 

In quasilinear regime particles oscillate about pancakes while drifting to- 
wards deep potential wells. Alignment of velocity and acceleration is disturbed, 
except in direction of bulk motion. In nonlinear regime, velocities are random 
and there is no alignment between velocities and acceleration. 
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2. Velocity Contrast 

Above discussion shows that velocity and acceleration are aligned before shell 
crossing, after this misalignment increases till velocities are randomized. There- 
fore misalignment of velocity and acceleration [ in direction and magnitude] is 
an indicator of nonlinearity. Velocity contrast quantifies this mismatch as 

_ (u-g)^ 

A comparison of velocity and density contrasts for CDM simulations is presented 
here. [ Shell crossing is necessary for formation of dynamically nonlinear objects, 
even if they form by merger of smaller virialised objects. This justifies use of 
velocity contrast for CDM spectrum.] This simulation used a PM code with 
128"^ particles in a 128"^ box, with Mpartide — 2 X W^^Mq, so results apply to a 
CDM universe smoothed at that [ mass] scale. 

Figure 1 shows contours of equal population on -Dgu ~ (1 + ^) plane for two 
epochs [ Z = 1,0]. Relation for spherical top hat model has also been plotted for 
comparison. Average relation between density and velocity contrast coincides 
with that expected from STH over a large range in both panels. At large values 
of velocity contrast [ -Dgu ~ 50], there is an upward shift of contours from Z = 1 
to Z = 0. At a given level of nonlinearity in dynamics, say at turnaround, 
density is highest for spherical and lowest for planar collapse. Therefore an 
upward shift in the relation implies a shift from predominantly planar collapse 
towards spherical collapse. Collapsing objects with different shapes contribute 
to the large dispersion. At low velocity contrast, dispersion arising from this 
effect is almost an order of magnitude [ See Bagla and Padmanabhan (1995)]. 

Comparison of velocity and density contrasts for a rich cluster is shown 
as contours of equal projected density and velocity contrast for a cube in fig. 2. 
There are distinct differences in two indicators, density contours always pick out 
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Figure 2. Contours of projected density and velocity contrast for a 
cube of size 14h~^ Mpc from a CDM simulation. 

regions that appear prominent to the eye. Velocity contrast does not always 
follow the same contours, at times there is no contour corresponding to a pan- 
cake. This happens if matter has not shell crossed and is only falling along a 
"pancake" . 

Velocity contrast indicates substructure in some density peaks, this results 
from recent merger of smaller virialised objects or accretion of matter with little 
nonlinear dynamics. Some dense objects are not visible in contours of velocity 
contrast. In such objects, bulk motion towards neighboouring density peaks 
dominates over random velocities within. Highest contour levels for the two 
indicators do not always coincide, which could arise from presence of infalling 
material that has not undergone sufficient mixing. 

Virialised structures, like galaxies can form only after shell crossing and 
subsequent mixing of matter. Therefore only regions with sufficiently nonlinear 
dynamics can host luminous structures. Clustering properties of these regions 
and underlying mass distribution can be used to define an effective bias : 



where ^ is the correlation function. In fig. 3 we have plotted the effective bias 
parameter as a function of scale for four threshold values of velocity contrast. 
These figures show that bias is larger for regions with greater nonlinearity in 
dynamics, and it approaches unity at large scales. 

Ability of -Dgu to delineate regions with different levels of dynamical non- 
linearity can be used for studying environmental effects in galaxy formation. 
Assuming that different types of galaxies form in regions of different levels of 
nonlinearity, we can compute relative bias for these. Figure 3, shows bias pa- 
rameter for three bins of velocity contrast. Bias parameter for mass in regions 
with .1 < Dbfgu < 1 is smaller than unity. This however could be an artifact of 
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Figure 3. Right panel shows bias for regions with -Dgu > -l, 1, 10 & 
20 respectively. Left panel shows bias for regions in intermediate bins 
.1 — 1, 1 — 10 & 10 — 20. Bias is higher for more nonlinear regions in 
both panels. 

granularity in simulations. Bias parameter for other bins is greater than unity, 
shows strong variation at small scales and approaches unity at large scales. 

3. Discussion 

Velocity contrast can delineate regions with different levels of nonlinearity in 
dynamics, this can provide new insight into nonlinear dynamics and clustering. 
As dynamics of coUisional matter is modulated by that of coUisionless dark 
matter, we can use simple modelling to understand some environmental effects. 
Velocity contrast can also be used as a preliminary filter for computationaly 
intensive schemes like DENMAX. ( Bertschinger and Gelb f991). 

A comparison of regions selected using velocity contrast with those where 
luminous structures form in full hydrodynamical simulations can be used for 
studying infiuence of dynamics of dark matter on baryonic structures and relative 
importance of astrophysical processes at various scales. 
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